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Abstract 

We examine highly super-Eddington black-hole models for SS 433, based on two-dimensional hydrody- 
namical calculations coupled with radiation transport. The super-Eddington accretion flow with a small 
viscosity parameter, a = 10 -3 , results in a geometrically and optically thick disk with a large opening 
angle of ~ 60° to the equatorial plane and a very rarefied, hot, and optically thin high-velocity jets region 
, around the disk. The thick accretion flow consists of two different zones: an inner advection-dominatcd 

~. ' zone and an outer convection-dominated zone. The high-velocity region around the disk is divided into 

, two characteristic regions, a very rarefied funnel region along the rotational axis and a moderately rarefied 

|"T i ■ high-velocity region outside of the disk. The temperatures of ~ 10 7 K and the densities of ~ 10~ 7 g cm~ 3 

in the upper disk vary sharply to ~ 10 8 K and 10~ 8 g cm -3 , respectively, across the disk boundary be- 
tween the disk and the high-velocity region. The X-ray emission of iron lines would be generated only in a 
confined region between the funnel wall and the photosphcric disk boundary, where flows are accelerated 
to rclativistic velocities of ~ 0.2 c due to the dominant radiation-pressure force. The results are discussed 
regarding the collimation angle of the jets, the large mass-outflow rate obserevd in SS 433, and the ADAFs 
\ and the CDAFs models. 
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(— i 1. Introduction 

SS 433 is a very unusual and puzzling X-ray source, which exhibits remarkable observational features, such as two 
oppositely directed rclativistic jets moving with a velocity of 0.26 c, its expected too-high energy, and a processing 
motion of the jets over a period of f62.5 d. Although a number of papers concerning observational and theoretical 
data on SS 433 have been published (e.g., Margon f984; Cherepashchuk f993), the true nature of SS 433 is still not 
clear. Before the discovery of SS 433, Shakura and Sunyaev (f973) already discussed the observational appearance 
of a super-Eddington accretion disk around black holes. Regarding these points, many studies have been made 
concerning the relativistic jets of SS 433 (Katz 1980; Meier f 982; Lipunov, Shakura 1982) and its collimated ejection 
(Sikora, Wilson 1981; Bcgelman, Recs 1984). Thus, the super-Eddington accretion disks were generally expected to 
possess vortex funnels and radiation-pressure driven jets from geometrically thick disks (Lynden-Bcll f978; Fukuc 
f982; Calvani, Nobili f983 ). The radiative acceleration of particles near to the disk was also studied by several 
authors (Bisnovatyi-Kogan, Blinnikov f977; feke 1980, 1989; Fukue f996; Tajima, Fukue 1998). A two-dimensional 
hydrodynamical calculation of a super-Eddington accretion disk around a black hole was first examined by Eggum et 
al. (1985, 1988), and was discussed regarding SS 433. They showed a relativistic jet formation just outside a conical 
photosphere of the accretion disk and a collimation angle of ~ 30° of the jet, where the dominant radiation-pressure 
force accelerates the jet to about 1/3 c. However, the jet's mass flux is only 0.4% of the input accretion rate and is 
too small for SS 433. 

Kotani et al. (1996) strongly suggested that SS 433 is a binary system under a highly super-Eddington regime 
of accretion, where the mass-outflow rate of SS 433 jets is restricted to be > 5 x 10 _7 M Q yr^ 1 , combining spectral 
X-ray observations of iron emission lines with hydrodynamical modeling. Brinkmann and Kawai (2000) have modeled 
the two-dimensional hydrodynamical outflow of the SS 433 jets, while developing their previous model of conically 
out-flowing jets, and discussed the relationship between the physical parameters of SS 433 jets and observationally 
accessible data. 

Blandford and Begelman (1999) investigated an adiabatic inflow-outflow solution of the accretion disk, where 
the super-Eddington flow leads to a powerful wind. King and Begelman (1999) discussed the super-Eddington mass 
transfer in binary systems and how to accrete in such super-Eddington accretion disks. Furthermore, recent theoretical 
studies and two-dimensional simulations of the black holes accretion flows ( Igumenshchcv, Abramowicz 1999; Stone 
et al. 1999; Narayan et al. 2000; Quataert, Gruzinov 2000; Stone, Pringle 2001; Abramowicz et al. 2001 ) showed 
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a new development in the advection-dominated accretion flows (ADAFs) and the convection-dominated accretion 
flows (CDAFs) which may be deeply related to the super-Eddington accretion flows. Okuda and Fujita (2000) have 
examined the super-Eddington accretion disk model of a neutron star for SS 433 by two-dimensional hydrodynamical 
calculations and discussed the characteristics of the super-Eddington flows. In this paper, following the neutron star 
model, we examine the super-Eddington black-hole models and discuss the results regarding SS 433 and also AVAFs 
and CDAFs. 

2. Model Equations 

A set of relevant equations consists of six partial differential equations for density, momentum, and thermal and 
radiation energy. These equations include the full viscous stress tensor, heating and cooling of the gas, and radiation 
transport. The radiation transport is treated in the gray, flux-limited diffusion approximation (Levermore, Pomraning 
1981). We use spherical polar coordinates (r,£,<p), where r is the radial distance, ( is the polar angle measured from 
the equatorial plane of the disk, and <p is the azimuthal angle. The gas flow is assumed to be axisymmetric with respect 
to Z-axis (d/dtp = 0) and the equatorial plane . In this coordinate system, the basic equations for mass, momentum, 
gas energy, and radiation energy are written in the following conservative form (Kley 1989): 

^+divM = o, (i) 



9{pv) ,. , . 
gt + div(pvv) = p 



(r - r g ) 2 



-^ + f r + divS r + -S rr , (2) 
or r 



— — h div(prwv) = -pi^tanC - ^ + div(rS^) + S^tanC + fa, (3) 

o ( otcosCv i 

dt — — + div(prcosCf v v) = div(rcosCSV), (4) 

d ps 

-qJ- + div(pev) = — p divv + <I> — A, and (5) 

^ + divFo + div(vE + v-P )=A- p ^ + a \ ■ F , (6) 

where p is the density, v = (v,w,v v ) are the three velocity components, G is the gravitational constant, M* is the 
central mass, p is the gas pressure, e is the specific internal energy of the gas, Eg is the radiation energy density per 
unit volume, and P is the radiative stress tensor. It should be noticed that the subscript "0" denotes the value in 
the comoving frame and that the equations are correct to the first order of v/c (Kato et al. 1998). We adopt the 
pseudo-Newtonian potential (Paczyhsky, Wiita 1980) in equation (2), where r g is the Schwarzschild radius. The force 
density / R = (f r ,f^) exerted by the radiation field is given by 

fn = P^F , (7) 

where k and a denote the absorption and scattering coefficients and Fq is the radiative flux in the comoving frame. 
S = {S r ,SQ,S v ) denote the viscous stress tensor. $ = (S V)u is the viscous dissipation rate per unit mass. 

The quantity A describes the cooling and heating of the gas, i.e., the energy exchange between the radiation field 
and the gas due to absorption and emission processes, 

k = pcK{S*-E ), (8) 

where S* is the source function and c is the speed of light. For this source function, we assume local thermal equilibrium 
5* = aT 4 , where T is the gas temperature and a is the radiation constant. For the equation of state, the gas pressure 
is given by the ideal gas law, p = R^pT/p,, where p is the mean molecular weight and Rq is the gas constant. The 
temperature T is proportional to the specific internal energy, e, by the relation p = (7 — l)ps = RopT/p, where 7 is 
the specific heat ratio. To close the system of equations, we use the flux-limited diffusion approximation (Levermore, 
Pomraning 1981) for the radiative flux: 

F = -^L-grad.E 0! (9) 

p[K + (T) 

and 



Pv = E -T Edd , (10) 
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Model 


M./M 


M,/M B 


M* ( g s- 1 ) 


a 


R* 


^max/ ^* 


BH-1 


10 


4 


8 x 10 iy 


lO" 3 


2r g 


220 


BH-2 


10 


4 


8 x 10 19 


0.1 


2r g 


220 



where A and T^dd are the flux-limiter and the Eddington Tensor, respectively, for which we use the approximate 
formulas given in Kley (1989). The formulas fulfill the correct limiting conditions in the optically thick diffusion limit 
and the optically thin streaming limit, respectively. 

For the kinematic viscosity, we adopt a modified version (Papaloizou, Stanley 1986; Kley, Lin 1996) of the 
standard a- model. The modified prescription for v is given by 

v = a c s mm.[H p ,H] , (11) 

where a is a dimcnsionless parameter , usually a=0. 001-1.0, c s the local sound speed, H the disk height, and H p = 
p/ | gradp | the pressure scale height on the equatorial plane. 

3. Numerical Methods 

The set of partial differential equations (l)-(6) is numerically solved by a finite-difference method under initial and 
boundary conditions. The numerical schemes used are basically the same as that described by Kley (1989) and Okuda 
ct al. (1997). The methods are based on an explicit-implicit finite difference scheme. N r grid points (=150) in the 
radial direction are spaced logarithmically, while grid points (=100) in the angular direction are equally spaced, 
but more refined near the equatorial plane, typically A£ = 7r/150 for ir/2 > ( > ir/6 and AC = 7r/300 for ir/6 > ( > 0. 

3.1. Model Parameters 

For the central star of SS 433, we assume a Schwarzschild black hole with mass M* = 10M© and examine the structure 
and dynamics of an accretion disk around the black hole and its surrounding atmosphere. From an observational 
constraint of the mass-outflow rate, M\ oss , of > 5 x 10 _7 M Q yr _1 (3 x 10 19 g s _1 ) in SS 433 (Kotani ct al. 1996), we 
adopt a mass accretion rate M* of 8 x 10 19 g s _1 , which corresponds to ~ 4Me, where Me is the Eddington critical 
accretion rate for the black hole, given by 

M E = --=, (12, 

where n e is the electron scattering opacity. 

For the viscosity parameter, a, we consider two cases, 10 -3 (BH-1) and 0.1 (BH-2). The inner-boundary radius, 
i?*, of the computational domain is taken to be 2r g and the outer-boundary radius, i? ma x, is selected so that the 
radiation pressure is comparable to the gas pressure at the outer boundary, where the Shakura-Sunyaev instability 
never occurs. The model parameters used are listed in table 1. 

3.2. Boundary and Initial Conditions 

Physical variables at the inner boundary, except for the velocities, are given by extrapolation of the variables near 
the boundary. However, for the velocities, we impose limited conditions that the radial velocities are always negative 
and the angular velocities are zero. If the radial velocity by the extrapolation is positive, it is set to be zero; that is, 
outflow at the inner boundary is prohibited. On the rotational axis and the equatorial plane, the meridional tangential 
velocity w is zero and all scalar variables must be symmetric relative to these axes. 

The outer boundary at r = i? ma x is divided into two parts. One is the disk boundary through which matter is entering 
from the outer disk. At the outer-disk boundary we assume a continuous inflow of matter with a constant accretion 
rate, M*. The other is the outer boundary region above the accretion disk. We impose free-floating conditions on this 
outer boundary (i.e., all gradients vanish) and allow for outflow of matter, whereas any inflow is prohibited here. We 
also assume the outer boundary region above the disk to be in the optically-thin limit, | -Z^b | — ^ cEq. This imposes a 
boundary condition on the radiation energy density, Eq. The initial configuration consists of a cold, dense, and optically 
thick disk and a hot, rarefied, and optically thin atmosphere around the disk. The initial disk is approximated by 
the Shakura-Sunyaev's standard model. The initial hot rarefied atmosphere around the disk is constructed to be 
approximately in hydrostatic equilibrium. 
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Fig. 1. Time evolution of the luminosity, L, for model BH-1, where is the Keplcrian orbital period at the inner boundary 2r g . 

4. Results 

4.1. Model BH-1 with a = 1(T 3 

The initial disk thickness, H/r, based on the standard disk model, is ~ 1 at the inner region and ~ 0.1 at the outer 
boundary. The ratio /3 of the gas pressure to the total pressure at the initial outer disk boundary is ~ 0.27 . We 
performed a time evolution of the disk until t = 4 x 10 Pd for model BH-1, where Pd is the Keplerian orbital period 
at the inner boundary. 

Figure 1 shows the time evolution of the total luminosity, L, for model BH-1. After an initial sharp rise of the 
luminosity, the luminosity curve descends gradually toward a steady state value. The luminosity L at the final phase 
is 1.6 x 10 39 erg s -1 and L/Le is ~ 1, which is far smaller than M*/Me(= 4), where Le is the Eddington luminosity. 

At the initial stage of the evolution, matter in the inner region of the disk is ejected strongly outward. Some of 
the ejected gas hit on the rotational axis and others propagate outward. The gas hitting on the axis leads to a high- 
temperature region along the axis and anisotropic radiation fields, in which equ-contour lines of the radiation energy 
density are concentric. This results in an outward radiation-pressure force, which dominates the gravitational force, 
and a high-velocity jet region is formed along the axis. At the beginning of the evolution, the high-velocity region is 
confined to a narrow region along the rotational axis, but gradually spreads outward from the axis. The spreading-jet 
gas interacts with the surrounding medium and finally settles down to a quasi-steady state. On the other hand, the 
initially optically thin atmosphere above the initial disk is occupied by the dense gas ejected from the inner disk with 
increasing time, and becomes optically thick at the final phase, forming a cone-like shape with a large opening angle 
of ~ 60° to the equatorial plane. We regard this optically thick and dense region as an accretion disk. The boundary 
between the high-velocity region and the disk is characterized by sharp gradients of density and temperature. 

In order to understand well the evolution of the disk and the high-velocity region, we animated the density and 
temperature evolutions. Figure 2 denotes nine snapshots of the density evolution, where the right-hand number in 
each snap shows the evolutionary time in units of Pj; also, a color legend for log p (g cm -3 ) is shown. The green and 
orange parts in the later snaps express the rarefied, hot, and optically thin high- velocity jets region, while the blue and 
dark blue parts denote the optically thick and dense disk region. In the high- velocity region, we distinguish the orange 
region (which we call the funnel region) from the green region, because the densities in the orange region decrease 
sharply towards the rotational axis and are as very low as 10 _18 -10 -26 g cm -3 , compared with 10 _8 -10 -18 g cm -3 in 
the green region. Therefore, we have three characteristic regions here: (A) the disk region (blue), (B) the high- velocity 
region (green), and (C) the very rarefied funnel region (orange). As discussed later, the centrifugal barrier and the 
funnel wall lie near the boundaries, between (A) and (B), and between (B) and (C) at the final phase, respectively. 

Figure 3 shows the density contours with velocity vectors (a) and the temperature contours (b) on the meridional 
plane at t = 3577Pd- In figure 3a, the density contours are denoted by the labels log p = —5, —6,-7, —8,-10, and —18. 
The velocity vectors show the maximum velocity, 0.4 c, at £ ~ 80° and r/R* ~ 200. The thick dashed and the thick 
dot-dashed lines in (a) show the funnel wall and the disk boundary between the disk and the high-velocity region, 
respectively. Here, we approximately define the disk boundary as an interface through which a temperature of ~ 10 7 
K in the disk jumps to ~ 10 8 K. The disk boundary lies roughly between the density labels —7 and —8. In the funnel 
region (C) ( C > 80°), the flow velocities are very large, but rather chaotic near the rotational axis. In the high- velocity 
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Fig. 2. Snapshots of the density evolution for model BH-1, where the right-hand number in each photo shows the evolutionary time 
in units of Pj; also, a color legend for log p (g cm -3 ) is shown. The first snapshot denotes the initial density configurations, which 
consist of the Shakura-Sunyaev disk and the surrounding rarefied atmosphere. The green and orange parts in the later snaps show 
the hot, rarefied, and optically thin high-velocity jets region, while the blue and dark blue parts denote the cold, dense, and optically 
thick disk region. The flow velocities at the final phase are ~ 0.1-0.2 c in the deep-green part and ~ 0.2-0.4 c in the light-green 
and orange parts. The observed X-ray emission of iron lines in SS 433 would be formed in the deep-green high-velocity region just 
outside the disk boundary. 
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Fig. 3. Velocity vectors and contours of the density, p (g s~ 3 ), in logarithmic scale (a) and bird's-eye views of the gas temperature 
T(K) (b), on the meridional plane at the evolutionary time t = 3577Pd for model BH-1. The density contours are denoted by the 
labels log p = —5, —6, —7, —8, —10, and —18, and the velocity vectors show the maximum velocity 0.4 c at f ~ 80° and r/ii* ~ 200. 
The thick dashed and the thick dot-dashed lines in (a) show the funnel wall and the disk boundary between the disk and the 
high-velocity region, respectively. 
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Fig. 4. Contours of the logarithmic radiation energy density log Eq in an arbitrary unit (a) and the logarithmic angular-velocity 
log Q normalized by the Kcplcrian angular velocity at the inner boundary with intervals of 0.25 (b) on the meridional plane at 
t = 3577Pd for model BH-1, where the thick dashed and the thick dot-dashed lines denote the funnel wall and the disk boundary 
between the disk and the high-velocity jets region, respectively. 

region (B) ( 80° > C ^ 60°), the flows are relativistic to be ~ 0.1-0.2 c, and the gas streams run radially. Figure 3b 
shows a bird's-eye view of the temperature, where the disk boundary is clearly recognized by a sharp wall with high 
temperatures. The temperatures range from 2 x 10 6 to ~ 10 7 K in the disk region, jump to ~ 10 8 K just across the 
disk boundary, and again distribute gradually between ~ 10 8 and 10 9 K in the high-velocity region. In the funnel 
region near the rotational axis, the temperatures are as very high as ~ 10 9 -10 11 K. On the mid-plane, the densities 
and the temperatures are not largely different from the initial ones, except for the innermost region of r/i?* < 10, 
where the flow is more rarefied and hotter than the initial flow. 

Figure 4 shows the contours of radiation energy density, Eq, in an arbitrary unit (a) and angular- velocity f2 (b) 
on the meridional plane at t = 3577Pd, where the thick dot-dashed line denotes the disk boundary. In figure 4a, Eq 
is almost independent on £ in the optically thin high-velocity region and weakly dependent on £ in the upper disk 
region, and so the radiation-pressure gradient forces act only radially in the high- velocity region (B). In figure 4b, the 
angular velocities, 0, in unit of the Kcplcrian angular velocity at the inner boundary are denoted by the contours of 
log f2. The angular velocities show two kinds of characteristic distributions. Near to the equatorial plane, the angular 
velocities are nearly Kcplcrian. However, in the upper disk and the high-velocity region, the equi-contours of Q are 
parallel to the Z-axis, and fl behaves approximately as K/x , where x is the distance in units of 2r g from the Z-axis 
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and K is a constant, although the constant K seems to change slightly on the disk boundary and the funnel wall. As 
a result, the specific angular momentum A a (= x 2 Q) is approximately conserved in these regions. From the numerical 
data, we find fl = 5.6/x 21 with our non-dimensional units in the high- velocity region. The funnel wall, a barrier where 
the effective potential due to the gravitational potential and the centrifugal one vanishes, is described by a surface 

where rf = {x 2 + z 2 ) 1 / 2 is given in units of 2r g (Molteni et al. 1996). From equation (13) together with the above O, 
the funnel wall is shown by the thick dashed line in figure 4b. The funnel region (C) and the high-velocity region (B) 
are roughly separated by the funnel wall. 

The acceleration for the relativistic flows depends on the distribution of the radiation pressure, P I (= /eEq), where 
/e is the Eddington factor. The flux limiter, A, almost equals to 1/3 in the optically thick disk, while in the optically 
thin high-velocity region, A is very small and its spatial variation is almost identical to that of the density, p. In 
regions (B) and (C), the radiation temperatures, T r = (Eo/a) 1 / 4 , depend almost only on r and are small compared 
with the gas temperature, T; also, the gas temperatures in the region (C) are much higher than that in the region (B). 
However, its too high temperature (~ 10 9 ~10 11 K) in region (C) may be unreliable because we did not take account 
of other physical processes, such as Compton processes and pair production-annihilation, which may be important at 
such high temperatures. On the other hand, we consider that the densities would not be drastically altered by the 
processes. 

The radiation pressures are very dominant everywhere, especially in the inner region of the disk. In the disk 
region (A), the centrifugal forces balance the gravitational forces near to the disk midplane, whereas, at the upper 
disk region, the balances due to the centrifugal, the gravitational, and the radiation-pressure gradient forces are 
maintained. In regions (B) and (C), the dominant radial radiation-pressure forces are one order of magnitude larger 
than the gravitational forces because the densities in these regions are very low. As a result, the radiation-pressure 
gradient forces accelerate the flows radially to relativistic velocities. The azimuthal velocities in the high-velocity 
region are in the range of 0.1-0.2 c, and the high-velocity gas is blown off spirally and relativistically. At the disk 
boundary between regions (A) and (B), the gravitational force balances roughly the radial component of centrifugal 
force, i.e., the boundary region corresponds to what is called a centrifugal barrier. 

Figure 5 shows magnified velocity fields and the density contours near the upper disk boundary at t = 3577Pd, where 
the maximum velocity is 0.22c and the contours of log p (g s -3 ) is shown by lines with labels —6,-7,-8, and —10. 
The disk boundary between regions (A) and (B) is between the density labels —7 and —8. As mentioned later, this 
high-velocity region above the disk boundary would correspond to a typical X-ray emitting region of iron lines with 
velocities of ~ 0.2 c. 

Figure 6 shows the flow features in the inner disk and the surrounding high-velocity region, where the velocities are 
indicated by unit vectors and the thick dot-dashed line shows the disk boundary. One of the remarkable features of 
luminous accretion disks is convective phenomena in the inner region of the disk (Eggum et al. 1985, 1988; Milsom, 
Taam 1997; Fujita, Okuda 1998). The convective motions are clearly found in this figure, and there appear more than 
a dozen of convective cells. 

Previously it has been shown that a luminous accretion disk with M* ~ Me is unstable against convection, and that 
this is induced by a large negative gradient of entropy vartical to the equatorial plane (Bisnovatyi-Kogan, Blinnikov 
1977; Fujita, Okuda 1998). The accretion flow is severely disrupted by the convective eddies. Recent theoretical 
studies and two-dimensional hydrodynamical simulations of radiatively inefficient black hole accretion flows (Stone et 
al. 1999; Igumcnshchev, Abramowicz 1999; Narayan et al. 2000; Quataert, Gruzinov 2000) showed that accretion flows 
with low viscosities are generally convection-dominated flows (ADAFs and CDAFs) and have characteristic self-similar 
solutions of the disk variables which are described by power law profiles with radius. This was also confirmed in recent 
magnetohydrodynamical simulations of the accretion flows around black holes (Stone, Pringle 2001). 

Model BH-1 belongs to a typical ADAFs or a CDAFs case because of its small value of the viscosity parameter, a. 
Figure 7 denotes the radial and angular profiles of the time-averaged flow variables near the equatorial plane for model 
BH-1: (a) density, p, in units of po(= 10 -8 g cm -3 ); (b) radial velocity, v, in units of Keplcrian velocity uk, at r = _R„; 
(c) total pressure, P t , in units of po v\ ; (d) rotational velocity, v v (solid line), in units of «k» and the Keplcrian one 
(dotted line); (e) radial mass-inflow rate, M- ln {r) (g s _1 ) (solid line), mass-outflow rate, M out (r) (dashed line), net 
mass accretion rate, M a (r)(= Mj n + M out ) (dotted line), in the whole region, and mass-outflow rate Mj et (r) (o) in the 
jet region; (f) convective luminosity, L c = J F c (r)dS (erg s _1 ); (g) entropy, S, at £ ~ 0° (solid line) and 42° (dotted 
line) in an arbitrary unit; and (h) angular profiles of mass-inflow rate, Mi n (£) (g s _1 ) (solid line), mass-outflow rate, 
Mout(C) (dashed line), and net mass-flow rate, M p (= M- ln + M out )(o), where F c in (f) is the convective energy flux 
given by Narayan et al. (2000) and its integral is taken within the disk region, and the ordinate scale in (h) shows 
Sign(M) log \M\. The time-averaged and angle-integrated mass inflow and outflow rates (Mj n and M out , respectively) 
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Fig. 5. Magnified velocity fields and density contours near the upper disk boundary at t — 3577Pd for model BH-1, where the 
maximum velocity is 0.22 c and the contours of log p (g s~ 3 ) is shown by lines with labels —6, —7, —8, and —10. This high-velocity 
region above the disk boundary is a typical X-ray emitting region of iron lines with velocities of ~ 0.2 c. 



in figure 7e are denned as 

Mi n (r) = Airr 2 I p min (v,0)cos (d(, 



and 



it/2 



M out (r) = 4irr 2 / p max (v, 0)cos C^C- 
Jo 



(14) 



(15) 



For these rates, the net mass accretion rate is M a = Mi n + M out . These mass- flux rates are actually regarded as those 
in the disk region, since Mj et is a few orders of magnitude smaller than M; n , M out , and M a in figure 7(e). Similarly, 
Min(C) and M ou t(C) are defined as 



Mj n (() = 4-7rcosC / r p min (w, 0)dr, 
Jr, 



and 



■^out(C) = 47rcos£ / r p max (w,0)dr. 



(16) 



(17) 



Abramowicz et al. (2001) suggests that low-viscosity accretion flows around black holes consist of two zones: an 
inner advection-dominated zone, in which the net mass- inflow rate, M- m , is small, and an outer convection-dominated 
zone, in which M; n increases with increasing radii. From figures 6 and 7, we can see that the accretion flow apparently 
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Fig. 6. Unit-velocity vectors and density contours of log p = —5, —6, —7, —8, —10, and —18 in the inner disk and the surrounding 
high-velocity jets region at t = 3577Pd for model BH-1, where the dot-dashed line denotes the disk boundary bewteen the disk and 
the high-velocity jets region. A transition between the inner advection-dominated zone and the outer convection-dominated zone 
occurs at r ~ 10i?*(= 20r g ). 

consists of two such zones. The transition radius, r tr , between the zones is ~ 10i?*(= 20r g ), which is larger than 10r g 
in Stone and Pringle (2001), but smaller than 30-50r g in Abramowicz et al. (2001). At r < r tr , convection is absent 
and the mass-inflow rate is very small. On the other hand, at r > r tr the flow is convectively turbulent and accretes 
slowly. 

From figure 7, we have rough approximations of p oc r, Pt r~ °' 85 , v v oc r -1 / 2 , and M ln oc r. The angle-averaged 
variables are averaged over twenty mesh points between < ( < 12°. However, it should be noticed that some variables, 
such as the entropy, may be considerably dependent on £, as is found in figure 7g, which shows S(( ~ 0) oc r~ 1,6 at 
r/R* > 12 but S(C = 42°) ~ constant at r/R* > 40. The entropies at C = and 42° sharply increase at r/R* < 12 
and 40, respectively, which lie almost in the disk surface, because of S oc T 3 / p in the radiation-dominated region. The 
entropy profiles within the disk show equi-contours vertical to the equator and considerable radial gradients near to 
the disk mid-planc, but radial equi-contours in the upper disk. 

Other radial profiles are compared with those of (1) p oc r -1 / 2 , P t oc r~ 3 / 2 , v oc r -1 / 2 , and M- m oc r for a self-similar 
solution (Igumcnshchev, Abramowicz 1999; Narayan et al. 2000), (2) p oc r -1 / 2 , P t oc r~ 3 / 2 , v = 0, v v oc r~ x / 2 , 
and Mi n =0 for a non-accretion convective envelope solution (Narayan et al. 2000; Quataert, Gruzinov 2000), and (3) 
p oc r , Pt oc r" 1 , v oc r~ 1 , and M- m oc r for hydrodynamical simulations with v = 10~ 2 p (Stone et al. 1999). The 
large differences of the profiles between us and them are the density profile with p oc r in model BH-1. However, it 
is natural that the very luminous disk like model BH-1 would have such density inversion profile, because the initial 
Shakura-Sunyaev disk has p oc r 3 / 2 and P t oc T 4 oc r~ 3 / 2 in the inner region of the disk, where the radiation-pressure and 
elctron scattering are dominant (Shakura, Sunyaev 1973). However, the self-similar solutions and other simulations 
are considered under a negligible radiation-pressure condition. 

In figure 6, we remark on the existence of an "accretion zone" just below the disk boundary, but above the top of 
the convective zones. This accretion zone is also found in Eggum et al. (1985, 1988). At r > rt r , roughly half the mass 
at any time at any radius will be flowing in and flowing out, respectively, and the mass-inflow rate, M- ln , balances 
the mass-outflow rate, M out . Near to the transition region, by way of the convective zones and the accretion zone, 
matter is accreted towards the equatorial plane. The accreting matter, which is carried to the transition region by 
convection, partly diverts into the high-velocity jets and partly flows into the inner advection-dominated zone. Thus, 
the mass-flow rate swallowed into the black hole is as very small as ~ 5 x 10 14 g s _1 , because the densities near the 
inner boundary arc as very small as ~ 10~ 14 -10 -26 g cm -3 , although the velocities arc as very large as —0.3c to —0.6c. 
On the other hand, the mass-loss rate of the jets at the outermost boundary is ~ 4 x 10 18 g s _1 , which is one order 
of magnitude smaller than the input accretion rate, M*, whereas the mass-outflow rate, M out , through the outermost 
disk boundary is as large as the input accretion rate. We notice that the radial mass-outflow rate in the jets region 
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Fig. 7. Radial and angular profiles of the time-averaged flow variables near to the equatorial plane for model BH-1: (a) density p , 
(b) radial velocity v, (c) total pressure Pt, (d) rotational velocity v v , (e) four radial mass-flux rates M(r), (f) convective luminosity 
L c , (g) entropy S, and (h) angular profiles of three mass-flux rates M(£) (see text in detail). 
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Fig. 8. Time evolution of the luminosity L for model BH-2. 

increases with radii, roughly as Mj e t oc r in figure 7e, analogously with other mass-flow rates M- m (r) and M out (r), and 
furthermore that the ^-direction net mass-outflow rate M p in figure 7h changes from its negative value to a positive 
value of ~ 10 18 g s _1 at £ ~ 56°, where the angle position corresponds to the boundary between the disk and the jets 
region. After crossing the boundary, M p decreases abruptly because the outflow gas from the disk surface is strongly 
bent toward the radial direction due to the dominant radial radiation-pressure force in the high-velocity region. 

This shows that a considerable disk wind is generated in the upper disk and is incorporated into the high-velocity 
flow. Therefore, a part of the convective flow escapes from the whole system as the disk wind and the other may 
always remain in convective circulations through the disk. Thus far, we expect that the mass-loss rate of the jets at 
great large radii will become much larger than the 4 x 10 18 g s _1 calculated here. 

4.2. Model BH-2 with a = 0.1 

The shape of the initial disk is very similar to that of model BH-1, but the densities and the temperatures on the 
disk mid-plane in model BH-2 are more than one order of magnitude and by a few factors larger than those in model 
BH-1, respectively. Figure 8 shows the time evolution of the luminosity for model BH-2. The luminosity, 2.5 x 10 38 erg 
s _1 , at the final phase, t = 4365Pd, is one order of magnitude smaller than that in model BH-1, and L/L^ is ~ 0.16. 

This model shows very different time evolutions from model BH-1. Mass ejection from the disk surface first begins 
at the innermost region of the disk, and subsequently occurs at the outer part of the disk. Although the ejected gas 
is blown off partly through the outer boundary, it accumulates in the outer part of the disk with increasing time. 
Finally, the disk becomes fatter at the outer part than the initial disk. The obtained accretion disk is geometrically 
thick as H/r ~ 0.2, and has an opening angle of ~ 10° to the equatorial plane. 

Figure 9 shows the density contours with velocity vectors (a) and bird's-eye views of the temperature contours (b) 
on the meridional plane at the final phase. The thick dot-dashed line in figure 9a shows the disk boundary between 
the disk and the surounding atmosphere. In this figure we never find any rclativistic jets along the rotational axis 
as is found in model BH-1. The temperatures and densities at the atmospheric region around the disk are as hot 
as ~ 10 8 ~10 10 K and as low as ~ 10 -9 -10 -12 g cm' 3 , respectively, but the outflow velocities are not very large in 
model BH-2, where the maximum velocity is ~ 0.07 c near the rotational axis. Large gradients of the density and the 
temperature across the disk boundary are also found in figures 9a and b. 

Figure 10 shows contours of the radiation energy density Eq (a) and the angular velocity Q (b) on the meridional 
plane at t = 4365Pd with the same units as in figure 4. The disk is nearly Keplerian throughout its whole shape, but 
the atmospheric region above the disk shows non-Kcplerian angular velocities. There also appear convective cells at 
15 t/R* ^ 50. The densities and temperatures at r/R* < 10 on the disk mid-plane become more rarefied and much 
hotter than those in the initial disk, respectively, and a spherical high-temperature region is formed around here. The 
densities near the inner boundary are much larger, ~ 10 -9 - 10~ 10 g cm~ 3 , than those in model BH-1. The resultant 
mass-inflow rate swallowed into the black hole is much larger, ~ —1.4 x 10 16 g s _1 , than the ~ —5 x 10 14 g s _1 in 
model BH-1, but is negligibly small compared with the input accretion rate. 

The disk features obtained here are interpreted, rather, under a category of the standard disk model, though the 
disk is geometrically thick. The result of model BH-2 describes that the concept of rclativistic jets expected for the 
super-Eddington accretion disks around black holes must be modified under the a-model of viscosity, i.e., for the case 
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Fig. 9. Velocity vectors and contours of the density p (g s~ 3 ) in logarithmic scale (a) and bird's-eye views of the gas temperature 
T (K) (b) on the meridional plane at t = 4365Pd for model BH-2, where the thick dot-dashed line in (a) shows the disk boundary 
between the disk and the surrounding atmosphere, and the flow vectors are denoted in the same units as in figure 3a. 
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Fig. 10. Contours of the radiation energy density (a) and the angular-velocity Q (b) on the meidional plane at t = 4365Pd for 
model BH-2 with the same units as in figure 4. 

of strong kinematic viscosity, there exists no rclativistic and well-collimatcd outflow jet, even in the super-Eddington 
accretion disks. Though the densities and temperatures at the atmosphere around the disk are in the range of those 
required for the X-ray jets of SS 433, the maximum outflow velocity is at most ~ 0.07c and the flows are not well 
collimatcd anywhere. Therefore, we can not expect X-ray emission lines with a definite Dopplcr shift of 0.26 c observed 
in SS 433, and conclude that model BH-2 is unfavorable for SS 433. 

4-3. X-ray Emitting Region of Iron Lines 

For a Maxwcllian distribution of the electron velocities, the power emitted per unit volume due to excitations of 
level n' of ion Z in the ground state n by electron collision is given by 

% = 1.9 X W- 16 T-V 2 n ( e- AE ' kT N c N z erg cm" 3 s"\ (18) 
av \ in / 

where AE, In, k, N c , Ni, and Q are the excitation energy between the n and n' levels, the ionization potential of 
hydrogen, the Boltzmann constant, the number density of electrons, the number density of ion Z, and the effective 
average of collision strength O, respectively (Blumcnthal, Tucker 1974). The red- and blue-shifted highly ionized 
iron lines, FeXXV K Q and Fe XXVI K Q , are clearly seen in all of the observed X-ray spectra of SS 433 (Kotani et 
al. 1996). Typically, focusing on the Fe XXV K Q line, we calculated the total power emitted by the iron lines through 
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Fig. 11. Distribution of mesh points contributed mostly to the total X-ray emission of iron lines at t = 3577Pd for model BH-1, 
where the disk boundary (thick dot-dashed line) and the temperature contours of log T = 7, 8, 9, and 11 (lines) are also plotted. 
The diamonds, triangles, filled triangles, and filled squares denote the mesh points with different flow velocities of ~ 0.10 c, 0.13 
c, 0.17 c, and 0.20 c, respectively. Only mesh points with emissivity of more than a hundredth of the maximum emissivity in the 
optically thin region are plotted here. 

the optically thin region. 

Figure 11 shows the mesh points which contributed most to the total X-ray power of the iron lines at t = 3577P c i 
for model BH-1, where only mesh points with emissivity of more than a hundredth of the maximum emissivity in the 
optically thin region are indicated by the triangles, squares, filled triangles, and filled squares, which show different flow 
velocities of ~ 0.10 c, 0.13 c, 0.17 c, and 0.20 c at the mesh points, respectively. The influences of photoionization and 
photoexcitation by the intense radiation field in the innermost region were not taken account of here. If these effects 
are considerable, the results in figure 11 may be invalid. In this respect, we mention that Brinkmann and Kawai (2000) 
showed the effect of photoionization on the iron-line flux ratio to be relatively small, using their jets model of SS 433 
and assuming that the emission at the base of the jets is given by a black body, where the temperatures and densities 
(4 x 10 s K and 10 _9 -10 -10 g cm -3 , respectively) used in their model are of the same orders as those in the marked 
region in figure 11. 

We notice that the emitting region of the iron lines in figure 11 is confined into a narrow region at 63° < £ < 68° just 
outside the disk boundary, and that the high-velocity gas in the funnel region would never contribute to the X-ray 
line emission, because the densities in the region are more than ten orders of magnitude smaller than those in the 
narrow region, regardless of the temperatures in the funnel region. The accretion disk, itself, docs not contribute to 
the X-ray line emission, because the temperatures are as low as ~ 10 6 -10 7 K and its emissivity decreases exponentially 
as e - AE / kT in spite of the high densities in the disk. Therefore, we consider that the unique velocity, 0.26 c, of SS 
433 would be attributed to the effective velocity in the above emitting region of the iron lines. 

5. Discussion 

5.1. Comparison with Other Results 

Our present results are compared with those of a super-Eddington black hole model by Eggum et al. (1985, 1988), 
who considered a black hole of M* = 3M Q and = 4Mg and used the flux-limited diffusion approximation for the 
radiation transport. The black hole mass is by a factor 3 smaller than ours, but the input accretion rate is the same 
as ours. They used a constant kinematic viscosity, v = 1.5 x 10 15 cm 2 s _1 , instead of the usual a-modcl for viscosity. 

There exist some qualitatively important and common results between our model BH-1 and their model, such as 
the relativistic funnel jets, the geometrically thick disk, the accretion zone, and the remarkable convective phenomena 
in the disk. In their results, the flows between the disk and the rotational axis are accelerated to relativistic velocities 
due to the radiation-pressure force. The collimation angle of the jets to the rotational axis is also as large as ~ 30° in 
model BH-1. The mass-outflow rate obtained by Eggum et al. (1985, 1988) is three orders of magnitude smaller than 
the input accretion rate, AT*, and accordingly the estimated Mi oss of SS 433. On the other hand, the mass-outflow 
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rate from the high-velocity jets in model BH-1 is much larger than that in Eggum et al. although it is one order of 
magnitude smaller than the estimated value in SS 433. Most of the accreting matter in Eggum et al. (1985, 1988) is 
swallowed into the central black hole through the inner boundary, but the mass-flow rate into the black hole in our 
model is negligibly small compared with the input accretion rate. The differences in the numerical results may be 
attributed to the high kinematic-viscosity, v, of 1.5 x 10 15 cm 2 s _1 in Eggum et al. (1985, 1988), which is more than 
two orders of magnitude larger than those at r/R* < 20 and at r/i?„ > 50, respectively, on the disk mid-plane in model 
BH-1. 

5.2. Convection-Dominated Accretion Flows with Small a 

As an alternative model to ADAFs models, since Narayan and Yi (1994), Blandford and Begelman (1999) proposed 
advection-dominatcd inflow-outflow solutions (ADIOS) which include a powerful wind that carries away mass, angular 
momentum, and energy from the accreting gas. A typical case with appropriate parameters shows that only a small 
fraction of the mass supplied will reach the black hole, accompanying a disk wind. In relation to these ADAFs, Narayan 
et al. (2000), Quataert and Gruzinov(2000), and Abramowicz et al. (2001) developed the convection-dominated ac- 
cretion flows (CDAFs) or a 'convective envelope solution'. The CDAFs models assume that the convective energy 
transport is primarily directed radially outward through the disk, but with no outflows, whereas in the ADIOS models 
the convective energy transport has a large ^-direction component toward the disk surface at high latitudes, which 
will drive strong outflows. 

The initial Shakura-Sunyaev disk for model BH-1 is originally unstable to convection, since the entropy, S ozT 3 /p, in 
the radiation-pressure dominant inner disk falls abruptly with increasing z. Convection tends to establish an isentropic 
structure along the Z-axis, which results in equi-contours of entropy vertical to the equator near to the disk mid-plane 
in model BH-1. The convective flows accompanied by a large mass-outfow from the disk surface may be relevant 
to the ADIOS models. On the other hand, large negative radial gradients of the entropy near to the mid-plane in 
model BH-1 also induce large convective energy transport directed radially outward through the disk, as can be seen 
in figures 7f and g. As a result, the accretion flows in model BH-1 may be dominated by a ADIOS model and also a 
CDAFs model. 

The convective envelope solution, which is special one of the CDAFs, was addressed by Narayan et al. (2000) and 
Quataert and Grunizov (2000). The convective envelope solution expresses that v ~ and the viscous dissipation 
rate, Q + , is ~ because the net shear stress via viscosity and via convection vanishes. This leads to no advection 
entropy, and thus the divergence of the convective energy flux vanishes from the energy equation, that is, F c {r) oc r~ 2 . 
The relation F c oc r~ 2 is also derived from the definition F c = —3a(c 2 /flK)pTdS/dr (Narayan et al. 2000) and using 
the radial profiles of p oc r and Pt oc r~°- 85 in figure 7, where c s is the sound speed. Figure 7f shows roughly L c ~ 
constant in the outer region; accordingly, F c oc r~ 2 if F c is independent of £. Even for net zero stress, the negative 
radial gradient of the entropy would be maintained if the high-entropy gas is lost in a considerable disk wind, as is 
found in model BH-1. Accordingly, model BH-1 may possess a circumstance near the convective envelope solution, 
although the mean radial velocity is not exactly zero. In this respect, further knowledge of the actual mechanism of 
the angular momentum transport by convection is required. 

Abramowicz and Igumenshchev (2001) suggest that CDAFs have a significant outward energy flux carried by 
convection, with a luminosity of L c = e c Mc 2 , where the efficiency, e c , is ~ 3 x 10~ 3 -T0~ 2 , independently of the accretion 
rate, and the radiative output comes mostly from the convective part. In model BH-1, the convective luminosity, L c , 
near the outer boundary is ~ 5 x 10 37 erg s _1 , which is far smaller than the radiative luminosity, L(~ 1.6 x 10 39 at t 
= 3577 Pd in model BH-1); this corresponds to e c ~ 10~ 3 , which agrees with that by Abramowicz and Igumenshchev 
(2001). However, the luminosity, L, comes mostly from the innermost hot region at r ~ r tr through the optically 
thin high- velocity region, instead of the convective-dominated disk. Indeed, M a GM/r tr gives 2 x 10 39 erg s -1 , which 
corresponds well to L, roughly using M a (r tI ) ~ 10 20 g s _1 and r tr ~ 20r g from figure 7. 

Finally, we emphasize that most of the simulations and theoretical models so far have been done on thin accretion 
disks without radiation, and that it is not still clear whether these results are valid for the much thicker and radiation- 
pressure dominant disks. Further theoretical studies on the radiation- pressure dominant ADAFs and CDAFs are 
required. 

5.3. Relativistic Jets, Collimation, and Mass- Outflow Rate 

Observations of SS 433 show many emission lines of heavy elements, such as Fe and Ni, which denote the unique 
red and blue Dopplcr shifts of 0.26 c in the X-ray region. The collimation angle of SS 433 jets seems to be as small as 
0.1 radian (several degrees) (Margon 1984). Any models for SS 433 must explain the characteristics of these emission 
lines. The unique velocity, 0.26 c, may be reasonably explained in terms of the relativistic velocities at a confined 
region just outside of the upper disk boundary, which lies between the funnel wall and the disk, as is found in model 
BH-1. The velocity 0.26 c is not an inherent velocity of SS 433 like an orbital velocity in the binary system, but is a 
result of the super-Eddington accretion disk. 
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The collimation angle, ~ 30°, of the high- velocity jets in model BH-1 is rather large compared with the value of ~ 
0.1 radian expected for SS 433 jets. If the calculated high- velocity jets should be truly confined to a small angle of ~ 
0.1 radian from the rotating axis, we need a collimation mechanism operating outside the present computation grids, 
or some mechanism, such as magnctohydrodynamical collimation, which is not considered here. This problem remains 
open, although there exist some discussions that the jets of SS 433 may not be well collimated, as is inferred from the 
line profiles, and that the 0.1 radian jet may not be consistent with the X-ray observations of SS 433 (Eggum et al. 
1988). 

Hydrodynamical modeling of the jets combined with recent X-ray observations reveals temperatures of 6-8 x 10 s K, 
particle densities of between 5 x 10 11 ~5 x 10 13 cm~ 3 at the base of the jets, and the length of the X-ray jets being > 10 10 
cm (Brinkmann ct al. 1991; Brinkmann 1993; Kotani et al. 1996). These temperatures and densities required for the 
X-ray jets qualitatively agree with those at the high-velocity region between the funnel wall and the disk boundary 
in model BH-1. The X-ray emitting region observed in SS 433 would be in the range of r = 10 10 -10 12 cm (Kotani et 
al. 1996), which far exceeds the present computational domain. The high-velocity jets region in model BH-1 would 
probably extend up to r ~ 10 12 cm. Additionally, from the characteristics of the results given in figure 3, we expect 
that such an extended high-velocity region would have temperatures and densities of ~ 10 s K and ~ 10~ 12 -10~ 10 g 
cm -3 . Although the mass-loss rate (4 x 10 18 g s _1 ) of the jets in model BH-1 is one order of magnitude smaller than 
the mass-loss rate ( > 3 x 10 19 g s _1 ) estimated in SS 433, we speculate that the actual calculated mass-loss rate at 
large radii may become comparable to the observed one, because the jets mass-flux would increase due to the wind 
gas from the disk at large radii. In such a case, we may have a problem that the jets gas may slow down below 0.26 
c if the enhanced region of the jets mass-flux is too far from the present computational domain. Alternatively, the 
super-Eddington models with a much larger input accretion rate, M*, than the present value, 4Mb, may be responsible 
for a mass-loss rate comparable to ~ 10 20 g s _1 . 

5.4. Central Object of SS 433 

Whether the central object of SS 433 is a black hole or a neutron star remains unsolved even at present, because 
we have not yet obtained any decisive evidence of a central source mass. From some observational and theoretical 
estimates of a massive mass and its too high energetics of SS 433, many astrophysicists seem to favor a black-hole 
hypothesis of SS 433 (Leibowitz et al. 1984; Fabrika, Bychkova 1990 ; Cherepashchuk 1993; Fukue et al. 1998; Hirai, 
Fukue 2001), while some people suggest that the compact object is a neutron star, from the viewpoint of a theoretical 
model and observations of the He 11 line (Begelman, Rees 1984; D'Odorico et al. 1991; Zwitter et al. 1993). 

In this respect, our present results are also compared with those of a neutron-star model by Okuda and Fujita (2000) 
with M* = 1.4M , i?* = 10 6 cm, and a mass accretion rate (M*) of 10 20 g s -1 , which corresponds to ~ IOOMe for 
a neutron star. The viscosity parameter, a = 10~ 3 , used is the same as that in model BH-1. The characteristic 
features in model BH-1 are also found in this neutron-star model, where there appear a high- velocity jet region with a 
collimation angle of ~ 10° to the rotational axis, the existence of a confined X-ray emitting region of iron lines in the 
high-velocity region, a geometrical thick disk with an opening angle of 80° to the equatorial plane, a mass-loss rate 
far less than 10 20 gs" 1 , and convective motions in the inner disk. Regarding the collimation angle, the neutron-star 
model with a smaller angle ~ 10° may be favorable for SS 433. However, except for this point, there is no intrinsic or 
large quantitative difference between the neutron-star model and the black-hole model; that is, we could not find any 
decisive proof for these objects as a candidate of SS 433. 

5.5. Assumption of the Flux-Limited Diffusion Approximation 

The radiation transport was treated here as an approximation of flux-limited diffusion. The flux limiter, A, and the 
Eddington factor, /e, used in the approximation are given from the empirical formulas fitted to some stellar atmosphere 
models (Klcy 1989). Recently, Turner and Stone (2001) addressed that the flux-limited diffusion approximation (FLD) 
is less accurate when the flux has a component perpendicular to the gradient in radiation energy density, and in optically 
thin regions when the radiation field depends strongly on the angle. The calculations considered in our case typically 
result in two optical regions, which are the optically thick disk region and the optically thin high-velocity region. The 
flux limiter, A, is almost 1/3 in the disk and, contrarily, very small everywhere in the high- velocity region. In an 
optically thick disk, FLD is sufficient. On the other hand, in the optically thin high-velocity region, the radiation 
fields are highly anisotropic and FLD may be less accurate. The radiation fields in the high-velocity region originate 
in radiation from the central hot region and the outer disk surface. However, the contours of radiation energy density 
are almost concentric around the central object. Actually, the radial radiation flux, which is attributed to the central 
hot region, is one order of manitude larger than the radiation flux normal to the disk surface. As a result, we consider 
that FLD in our code would not result in any serious influences on the whole dynamics of the flows and the total 
luminosity; still, this problem should be checked with a more accurate method in future improvements of the numerical 
code. 
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We examined highly super-Eddington black-hole models for SS 433, based on two-dimensional hydrodynamical 
calculations coupled with radiation transport. A model with a small viscosity parameter, a = 10 -3 , shows that a 
geometrically and optically thick convection-dominated disk with a large opening angle of ~ 60° to the equatorial 
plane and rarefied, very hot, and optically thin high-velocity jets region around the disk are formed. The thick 
accretion flow near to the equatorial plane consists of two different zones: an inner advection-dominated zone, in 
which the net mass-inflow rate, M- m , is very small, and an outer convection-dominated zone, in which M- lri increases 
with increasing radii. The high-velocity region along the rotating axis is divided into two characteristic regions by the 
funnel wall, a barrier where the effective potential due to the gravitational potential and the centrifugal one vanishes. 
A confined region in the high-velocity jet region just outside of the photosphcric disk boundary would be responsible 
for the observed X-ray iron emission lines with a Dopplcr shift of 0.26 c. However, from this model, we can not obtain 
a small collimation degree of ~ 0.1 radian of the jets and a sufficient mass-outflow rate of the jets comparable to 
~ 10 20 g s^ 1 , as is expected for SS 433. These problems still remain open. On the other hand, from a model with a 
large a = 0.1, we find a geometrically and optically thick quasi-Keplerian disk with an opening angle of ~ 10° to the 
eauatorial plane, whereas the disk is far thinner than the disk with a = 10~ 3 . This model may be unfavorable to SS 
433 because we never find relativistic jets with ~ 0.2 c here. Further investigations of the super-Eddington models 
with other model parameters and over a wider range of the computational domain ( r ~ 10 12 cm) may give some key 
to the open problems. 

The author gratefully thanks Dr. D. Molteni for his useful discussions of the numerical method and the results, 
and also the referee for many useful comments concerning the advection-dominated and the convection-dominated 
accretion flows. Numerical computations were carried out at the Information Processing Center of Hokkaido University 
of Education. 
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